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AN INVESTIGATION OF WING BUFFETING RESPONSE 
AT SUBSONIC AND TRANSONIC SPEEDS: 

PHASE I F-111A FLIGHT TEST DATA ANALYSIS 

VOLUME I - SUMMARY OF TECHNICAL APPROACH, RESULTS AND CONCLUSIONS 

by 

David B. Benepe, Atlee M. Cunningham, Jr., 
and W. David Dunmyer 

SUMMARY 

The structural response to terodynamic buffet during moderate 
to high-g maneuvers at subsonic and transonic speeds was investi- 
gated. The measurements which consisted of shear, bending moment 
and torque at four wing span stations, vertical accelerations at 
the wing tips, center of gravity and pilot's seat and lateral 
accelerations at the center of gravity and pilot's seat had been 
previously obtained during the Loads Demonstration flight program 
on a variable sweep fighter- bomber aircraft. 

Existing flight data for one wing sweep were extracted from 
magnetic tape records and subjected to statistical analyses to 
determine the power spectra and root-mean-square values of re- 
sponse for each of the measurements at several angles of attack 
for each of seven maneuvers. The frequency content of the vari- 
ous responses is correlated with results of ground vibration 
tests to identify the response modes. The rms values of response 
are plotted against angle of attack to show the variations of 
intensity of response during the maneuvers. 



The investigation showed that the structural response to 
•i ffet is very complex. Almost all of the natural vibration 
modwS of the aircraft, both symmetric and antisymmetric, are ex- 
cited during buffet encountered in a moderate to high-g maneuver. 
Some of the sensors show pronounced changes in the relative modal 
contributions to the total response as penetration beyond buffet 
onset increases. The fluctuating shear and bending moment loads 
on the wing are small in terms of design loads except near the 
wing tip. The wing structural response in torsion is larger than 
anticipated on the basis of previously published buffet studies 
and amounts to between 1/3 and 1/2 of the rms values of wing 
bending response at high angles of attack. For the particular 
aircraft geometry and structure examined, there is some evidence 
of bending- tors ion coupling starting at angles of attack between 
9 and 12 degrees at Mach numbers of 0.80 and 0.87. 

The investigation is reported in three volumes. NASA 
CR-152109 presents a summary of tne technical approach, the results 
and conclusions drawn from the results. NASA CR-152110 presents 
plotted variations of Power Spectral Density (PSD) data with fre- 
quency for each structural response item for each data sample 
analyzed during the course of the investigation. NASA CR-152111 
presents Power Spectral Density (PSD) data in tabular form for the 
convenience of those who might wish to perform additional analysis. 
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SYMBOLS 


Note: Quantities are presented in the International System of 

Units (U.S. customary units in parenthesis). The work 
was performed using U.S. customary units. 

b wing span - m, (ft) 

B.M.ogs design value of wing bending moment. N-m. (in - lb) 
c.g.,C.G. "center of gravity" 
f frequency, hertz 

f Q spectral base frequency or analysis bandwidth, hertz 

wing vertical shear as measured by strain gages - 
N, (lb) 

g gravitational acceleration 

M Mach number 

M x Wing Bending Moment as measured by strain gages N-m, 

(in - lb) 

My Wing torsional moment - N-m, (in - lb) 

n m ax maximum maneuver load factor - g's 

S theoretical wing are (leading and trailing edges of 

swept panel extended to airplane centerline m , (ft^) 

T length of input frame in spectral analysis - seconds 

start time of interval for spectral analysis - seconds 

T2 stop time of interval for spectral analysis - seconds 

4 T time interval used for spectral analysis “ T2-T^,sec 

V DES design value of wing vertical shear, N, (lb) 

y lateral acceleration g's 

z vertical acceleration g's 
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SYMBOLS (Continued) 


a indicated angle of attack referenced to wing manufac- 

turing chord plane 


a max 

a nom 

a 

Aa 

0 

9 a 
" V max 

Xax 

*T 


maximum indicated angle of attack - deg. 

nominal angle of attack representing time interval AT 

indicated angle of attack at time Tj, deg 

increment in indicated angle of attack during time 
interval AT, deg 

indicated sideslip angle, deg 

rms value of acceleration fluctuations - g, rms 

maximum rms value of wing vertical shear fluctuations 
N, rms, (lb, ms) 

maximum rms value of wing bending moment fluctuations 
N-m, ms, (in - lb, rms) 

average ras value determined from power spectral 
analysis 
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ABBREVIATIONS 


i 

'i 

j 

Alt altitude 

Asy® antisymmetric 

B.M. bending moment 

t 

. *■ Cross -PSD, XPSD Cross power spectral density • 

dB decibel 

Dyn Press dynamic pressure 

FM frequency modulation - 

H z hertz 

t<or,hori horizontal 
in-lb, IN-J.B inch-pound 
inb*d inboard 

f 

L left 

lb, LB pound 

L/H left hand 

LWT left wing tip 

m meter 

N newton 

N-m,N-M newton-meter 

outbd outboard 

P.S. pilot seat 

PSD power spectral density 

R right 

right hand 

!> i 


R/H 



ABBREVIATIONS, (Continued) 


rms 

RWT 

Svni 

TOR 
w.s . 


root -mean-square 
right wing tip 
symroet ric 
torsion 

Wing Station for strain gage measurements 
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SECTION 1 


INTRODUCTION 

The phenomenon of aerodynamic buffet has been a challenge to 
aircraft design teams for many years. With the advent of truly 
high performance fighter aircraft which are capable of operating 
at high angles of attack, the intensity of buffet and the magni- 
tudes of the aircraft structural responses to buffet have become 
important design considerations. 

The state-of-the-art of buffet research is such that devel- 
opment of a valid method of predicting aircraft structural re- 
sponse to buffet appears feasible in the near future. A signifi- 
cant problem that exists is the dearth of published flight data 
to provide an adequate test of the validity of a prediction 
method. 

The investigation described in this report is an attempt to 
supply data of sufficient scope in terms of the number and types 
of flight measurements and in terms of the depth of analysis of 
the measurements to use for correlation with predicted response 
characteristics. A secondary, though no less important, objective 
of the investigation is to add to the aeronautical community's 
understanding of buffet phenomena which has advanced rapidly in 
recent years. 
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Considerable effort has been expe ded in research programs 
conducted by the National Aeronautics and Space Administration's 
Ames, Langley, and Flight Research Centers (Refs. 1-7), by the 
armed services (Refs. 8-14), and by airframe manufacturers (Refs. 
15-17). Each of these previous efforts has contributed something 
of significance to our understanding of buffet and of aircraft 
responses to buffet. 

Results of various flight test programs have shown that a 
most reliable indicator of maneuvering buffet onset is the abrupt 

change in rms response of a wing- tip accelerometer. It can also 
be inferred that wing-tip accelerometer response is a good indi- 
cator of the variation of buffet intensity with angle of attack. 
However, measurements need to be obtained of accelerations at 
several points on an aircraft and of the spanwise distributions 
of dynamic structural loads on the wing to obtain a true test 
of a prediction method. Additionally, the power spectra data 
are needed for each measurement over a range of frequencies 
which covers the important natural vibration modes of the air- 
craft structure. 

Recent developments in flight test instrumentation and data 
recording and data processing systems permit the use of power- 
ful techniques of random data analysis to study the buffet prob- 
lem. Power spectral density (PSD) plots can be readily constructed 
to obtain the needed variations of response with frequency and 
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calculate root -mean-square values (rms) of the accelerations and 
loads to obtain statistical measures of the variations of re- 
sponses with angle of attack. 

The duration of the flight maneuvers is short, therefore the 
data samples are usually not long enough to strictly satisfy the 
mathematical criteria for obtaining a high level of confidence in 
the spectral or nns estimates. Nevertheless, spectral estimates 
and nns values are much more meaningful quantities than the peak- 
to-peak or half -amplitude measurements prevalent in early studies 
of the buffet problem. 

It is appropriate at this point to describe the philosophy 
of this report. The main body of this report serves to describe, 
in detail, the techniques used in obtaining the flight data, the 
instrumentation, the data reduction and data analysis procedures 
and discuss and interpret selected results which illustrate the 
major findings of the investigation. Included in the appendix 
are time histories of pertinent data items in plotted fotm. 
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The flight program from which data have been extracted for 
use in this study was primarily for F-lllA flight loads demon- 
stration. The test aircraft was therefore instrumented with 
. numerous strain-gage sensors and accelerometers. Flight maneuvers 
included wind-up turns, pull-ups and roller coasters at various 
altitudes, Mach numbers, and target load factors. While it was 
not a particular goal of the flight program to investigate buf- 
feting, many of the flight test conditions were such that signif- 
icant levels of buffeting occurred. An opportunity thus existed 
to subject the extensive flight data to analyses specifically 
aimed at providing detailed buffeting response characteristics. 

The basic approach used in the study is described in the 
following paragraphs. 

The various records pertinent to the Flight Loads Demonstra- 
tion program were surveyed to identify candidate flight points 
for the investigation. About 90 combinations of wing sweep, Mach 
number, altitude and target load factor were selected for Initial 
investigation. Analog strip-chart playouts of thirty items of 
instrumentation were then made for each candidate flight point. 

The strip charts were used primarily to identify the particular 
flight maneuvers in which the wing responses were of sufficient 
magnitude to be meaningful in the study and to eliminate points 
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in which excessive wing spoiler activity mi^ht "contaminate" the 
response data. 

Since the investigation was aimed at providing flight data 
for comparison with predictions an additional criterion in the 
selection of flight points for detailed analysis was the exis- 
tence of corresponding wind tunnel data in terms of wing sweep 
and Mach number. The investigation finally was concentrated on 
seven flight maneuvers all of which were performed at one nominal 
wing sweep. Existing digital data for the selected flight points 
were then reviewed and plots were made of the variations of Mach 
number, angle of attack and load factor as functions of time dur- 
ing the particular maneuvers. Specific time intervals corre- 
sponding to nominal average angles of attack were selected for 
each maneuver and stochastic (statistical) analyses were per- 
formed on the outputs from seven accelerometers and twelve 
strain-gage sensors for each selected time interval. 

The major effort in the stochastic analysis was devoted to 
obtaining power spectral density (PSD) plots as a function of 
frequency and root -mean- square (rms) - rlues of the magnitudes of 
wing shear, bending moment, and torque at four wing span stations 
and accelerations at the wing tips, center of gravity and pilot's 
seat. These PSD and rms data are essential for correlation with 
prediction methods. In addition the PSD plots were used to 
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identify the frequencies and magnitudes of the major responses 
for each time interval. The particular vibration modes associ- 
ated with the major responses were then identified insofar as 
possible from the results of extensive ground vibration tests 
performed at an earlier time on F-111A aircraft (References 18 
and 19) . 

The variations of rms values of the various outputs (for 
specific ranges of frequencies) of instrumentation were plotted 
against angle of attack and comparisons made to observe trends 
with changes of altitude or Mach number. 
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SECTION 2 


AIRCRAFT DESCRIPTION 

The test aircraft was F-111A Number 13. A drawing showing 
the general features of the aircraft is presented as Figure 1. 
Detailed geometry associated with the aircraft and its components 
appears in Table 1. The aircraft has a variable sweep wing and 
a convention was adopted early in the development program that 
all aerodynamic coefficients would be referenced to geometric 
characteristics at a specific wing sweep, namely, Ajjg * 16 de- 
grees. The variations of some key geometric characteristics of 
the wing with wing leading-edge sweep angle are presented in 
Figure 2. 

Although the aircraft is fitted with a high lift system 
consisting of multisegment leading-edge slats and multisegment 
double-slotted trailing-edge flaps, these devices were in their 
retracted positions for all maneuvers analyzed in this study. 

TVto-segment upper surface spoilers on each wing are used at 
low wing sweeps in addition to differentially controlled all- 
movable horizontal tails to achieve roll control. 

The aircraft has a three-axis stability augmentation system 
which was operational on all maneuver b analyzed in this investi- 
gation. 
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Figure l. F-lllA THREE- VIEW 
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SECTION 3 


AIRCRAFT INSTRUMENTATION 

The instrumentation system installed ir the aircraft consis- 
ted of two 30 track and one 14 track FM analog magnetic tape re- 
corders and various transducers throughout the airplane. IRIG B 
time reference signals were recorded on each tape recorder to 
provide time correlation. The general locations of the accelero- 
meters pertinent to the buffet study are shown in Figure 3. The 
actual locations in terms of aircraft geometry references are 
listed in Table 2. 

The characteristics of the accelerometers most of which were 
commercially available units, are indicated in Table 3. The accu- 
racies quoted refer to the nominal flat frequency response up to 
the limit frequency quoted. No calibration data exist above the 
quoted limit of flat fre uency response, however, the natural 
resonant frequencies are well beyond 100 hertz for all of the 

accelerometers. There is no reason to suspect a significant 
deviation from the basic calibration factors quoted for 
frequencies up to 100 hertz. 


The locations of the strain gage sensors pertinent to the 
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buffet study are shown in Figure 4. Shear, bending moment and 
torque were measured at each of the four indicated wing stations 
on the right wing. The sensitivities of these measurements were 
governed by the fact that the wing loads were to be measured 
during maneuvers at load factors up to the maximum capability of 
the aircraft. As a consequence the signal-to-noise ratios for 
the present buffet studies were lower than is desirable. The 
calibration slopes for each channel of information are shown in 
Table 4. 

In several cases the frequency response upper limit for the 
measurements was set by the subchannel characteristics of the 
flight recording system. Table 5 lists the appropriate nominal 
limit of frequency response for each item of instrumentation 
based on the recorder subchannel arrangements for each flight 
selected for detailed analysis. 

Correlating items such as angle of attack, airspeed, Mach 
number, altitude, gross weight, and control surface position 
indication of the spoilers and horizontal tail surface deflec- 
tions were also recorded on the FW tapes. 

A special test nose boom was fitted to the aircraft to 
obrain angle of attack, sideslip angle, altitude and Mach number 
data . 
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SECTION 4 


BASIC DATA PROCESSING METHODS 

During the Loads Demonstration Flight Program, the FM analog 
magnetic tapes containing raw flight test data were processed by 
automated techniques. The real time data were first disp yed 
on strip chart recorders for instrumentation verification. Next, 
the data were digitized at sample rates of 1 to 20 samples per 
second under computer control. The specific sample ratio de- 
pended on user group requirements. The digitized data were 
then scaled, calibrated and output in computer listings and 
computer tapes for additional processing on an IBM System/360. 
Second generation computer runs were made to obtain corrected 
flight condition data such as gross weight, Mach number, alti- 
tude, dynamic pressure and fuel distribution at 1-second intervals. 

Microfilm records of the computer listings from the original 
flight program data reduction were used in the present program to 
make plots of angle of attack, normal load factor, Mach number 
and dynamic pressure as functions of flight time and to identify 
the gross weights and altitudes for the selected flight maneuvers. 
The Mach number, altitude and dynamic pressure data include cor- 
rections for position error. The angles of attack are indicated 

angles and do not include the effect of upwash at the nose boom. 

* 
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The following formula can be used to obtain an estimate of the 
true geometric angle of attack values if desired: 

= 0.318 + 0.931 a(degrees) 

This correction was not applied to the data presented in this 
report because the magnitude is not large in the range of angles 
of attack covered by the flight data and is within the uncertainty 
of the flight measurements taken under buffeting conditions with 
an aeroelastic aircraft. 

Time histories were made of about 30 items of instrumentation 
which were considered pertinent to the buffet study. Examples of 
each of the strip chart records are presented in Figure 5. Addi- 
tional time histories for the selected flight maneuvers are con- 
tained in the appendix. The particular groupings of items on the 
strip charts permit assessment of the time phasing of the build- 
up of the wing responses to buffet as angle of attack increases 
during the maneuver. The chart speed was nominally 20 mm/sec and 
the frequency response limit of the pen recorders is about 50 
hertz. For correlation purposes, a time reference trace with pen 
deflections each second appears in the lower border of the strip 
chart. The flight times associated with the particular maneuvers 
are annotated below the time trace. 

The Mach number, altitude and angle of attack traces repre- 
sent uncorrected indicated values. Occasionally these items 
show finite jumps from one bound of the record to the other. 
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Figure 5. TIME HISTORIES OF INSTRUMENTATION OUTPUT FOR WINDUP TURN 
FLIGHT 77, RUN S&C-R 
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Figure 5. Continued 
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Figure 5. Cont 


30 


>nrttttnttrr 

h 

<•) 



l tl 4 ntt'tr 




I* I' 


Figure 5. Concluded 



This situation occurs because the calibration is a anil ti -stepped 
function and th» measurements cross the boundaries of one or 
more ranges of values. 
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SECTION 5 


FLIGHT CONDITIONS SELECTED 
FOR DETAILED ANALYSIS 

The tune histories for the 90 candidate flight Maneuvers 
were studied thoroughly to deternine which Maneuvers Might pro- 
vide a set of flight buffet data that would serve as a good test 
of a prediction Method. Five basic criteria were used in the 
selection process: 

First, the wing responses should be large enough 
to Make a detailed analysis meaningful. 

Second, the nominal flight conditions should closely 
approximate Mach number and angle of attack conditions 
which had already been obtained in a wind tunnel test 
program. 

Third, effects of Mach number (flow field) and altitude 
(aeroelasticity) should be included, if possible. 

Fourth, the selected maneuvers should be free 
from effects of spoiler deflection, \f possible. 

Fifth, one of the selected points should closely 
approximate a flight point that had already been ob- 
tained in a NASA/FRC flight test program which used 
F-lllA Number 6 (Reference 5). 


33 



The seven maneuvers selected for detailed analysis are 
listed in Table 6. They consist of two wind-up turns, four pull 
ups and one roller coaster, all performed at a nominal wing sweep 
of 26 degrees. The Mach numbers are approximately 0.70, 0.80, 
and 0.87 for the three high altitude maneuvers (above 6000 
meters). At the two higher Mach numbers the selected maneuvers 
are for three altitudes, nominally 1500, 3700 and above 6000 
meters. The gross weights range from 266,000 II to 330,500 M. 

Variations of angle of attack, load factor, Mach number and 
dynamic pressure with flight time are presented in Figure 6 for 
each of the selected maneuvers. Inspection of Figures 6(a) 
through 6(g) shows that the wind-up turns are gradual maneuvers 
whereas the pull ups and roller coaster are rather abrupt maneu- 
vers. As a consequence, it was feasible to select more points 
within each of the wind-up turns for stochastic analysis than 
for the other maneuvers. 

Table 7 lists the segments of each maneuver selected for 
detailed analysis. In most cases the time duration of the rec- 
ords is one second, but some longer records were used particu- 
larly for the wind-up turn of Flight 77, S&C-R. The table also 
lists the indicated angle of attack at the start of each record 
(«l), at the end of each record (<> 2 ) «nd in a few cases the maxi- 
mum indicated angle of attack (*max) occurring during the record. 
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Figure 6* Concluded 
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Also presented is the increment in angle of attack during the 
record (Ao) which is relatively large for the pullup and '•oiler 
coaster Maneuvers . A nominal angle of attack ("non) has been 
assigned to each segment which is used later to plot trends in 
the variations of instrument responses with angle of attack. 

An adjunct to the point selection process was the dete ruina- 
tion of the buffet onset characteristics for each of the Maneu- 
vers from the tiae histories. The criterion used for bulfet 
onset was the first detectable change in level of response froa 
any of the wing instruaentation which was followed by an ever 
increasing level of response. The results of the buffet onset 
determinations are presented in Figure 7. It was possible to 
obtain buffet onset points for six of the seven maneuvers. The 
aircraft was slightly into buffet even at 1-g trim conditions 
for Flight 78, Run 4. Figure 7(a) shows that the variation of 
angle of attack for buffet onset with Mach number is saooth de- 
spite the significant differences in dynaaic pressure and pitch 
rates for the various Maneuvers. Figure 7(b) shows that the 
normal force coefficients for buffet onset show a slight trend 
of increasing with a decrease in altitude (increase in dynaaic 
pressure). This result is likely due to the fact that less 
horizontal tail deflection is required to produce the maneuvers 
and the loss in lift due to tail deflection is then smaller at 
the lower altitudes. 
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SECTION 6 


STOCHASTIC ANALYSIS TECHNIQUES 

The analysis techniques used in this study are compatible 
with American National Standard (ANS S2. 10-1971) recommended 
methods for analysis and presentation of shock and vibration 
data. A quick- look examination was performed on each time-his- 
tory measurement to determine the data classification, degree of 
stationer! ty, record length, and recoverability. 

Measurements 

Data reduction was performed on the following data: 

1. Shear, bending moment, and torsion at four wing 
stations (12 measurements). 

2. TWo wing-tip accelerometers (vertical) 

3. TWo C.G. vertical anti one C.G. lateral accelero- 
meters . 

4. Pilot seat vertical and lateral accelerometers. 

The analysis performed on these items consisted of PSD (includ- 
ing an average rms value), cross-PSD, and Probability Density 
and was distributed as follows: 

PSD (all items) 779 plots 

Cross-PSD (AM001 and AM002) 7 photographs 

Probability Density 126 photographs 
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In addition to the plots and photographs, a digital magnetic 
tape recording was made of spectral coefficients of all PSD's 
and Cross-PSD's. The magnetic tape format (shown in Fig. 8) for 
a single PSD consists of two adjacent records; the first record 
for identification and the second for raw data. 

Special-Purpose Processing 

A block diagram of the special-purpose stochastic equipment 
is shown in Fig. 9. The FM signal is discriminated to recover 
the analog signal. Band-pass filters at 3 Hz and 100 Hz (48 dB 
per octave) were used to reject unwanted frequencies and to 
minimize aliasing effects on the sampled data. The data is 
calibrated at this point and converted to non-dimensional quanti- 
ties. The l/D 100 Analyzer was used to compute the PSD's, cross- 
PSD's, and Probability Densities. The stochastic algorithms 
utilized by the T/D 100 to perform these functions are discussed 
below. 

After each PSD cross-PSD was plotted, the spectral coef- 
ficients were punched, in ASCII format, on perforated paper 
tape. This paper tape was then fed into an SEL-810A mini-com- 
puter for transcription to magnetic tape. This was done to in- 
crease storage density, for ease of handling, and for speed of 
subsequent playback. 
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Figure 8, STOCHASTIC DATA TAPE FORMAT 
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Auto-Spectral Density (PSD) 

The T/D 100 computes the PSD coefficients by first approxi- 
mating the complex Fourier transform of the input signal. The 
Fourier transform of the time-domain input function x(t) is given 
by: 



x(t)(cos 2wft - j sin 


2nft) dt 


( 1 ) 


where j » ^-1 . Since the time-domain input is sampled and 
quantitized in the analyzer, and only a finite number of samples 
are available, the finite transform is used, and separated into 
its real P(f) and imaginary Q(f) components can be written 

as follows: 

.112 

Px(f) « / x(t) cos 27rft dt (2) 

-T/2 

1/2 

Q T (f) - / x(t) sin 2jrft dt (3) 

•'-T/2 


where T is the length of the input frame, which is assumed to be 
centered about time t“0. 

Replacing the continuous input, x(t), with a set of 2N+1 
discrete samples at intervals of t 0 - and replacing the 
sinusoidal functions by corresponding values, the continuous 
integrals may be expressed as the sum of products: 
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(4) 


+N 

P(kf 0 ) * /J x(nt 0 ) cos [2 
n«-N L 



+N 

Q(kf c ) - x(nt 0 ) sin [2 kf 0 (nt 0 )l (5) 

n»-N L J 


where k is a series of 2N integers and f Q is the base frequency 
which is equal to 

The PSD coefficients £s(kf 0 )J are then computed from (4) and 
(5) by the equation: 


S(kf D ) = 


P(kf 0 ) 



Q(kfo) 


( 6 ) 


Cross -Spectral Density (XPSD) 

The T/D 100 computes the XPSD (S A g) by simultaneously com- 
puting the complex Fourier transform of channel A (P A + jQ A ) and 
channel B (Pg + jQg) and multiplying the channel A transform by 
the complex conjugate of the channel B transform as follows: 

S AB * < P A + iQ A H p B " JQb> 

" P A P B + QaQb + J(P B Q A - P A Q B ) (7) 

The spectrum is then separated into its real (S p . AB ) and imagi- 
nary (Sq. ab ) components as follows: 

S P* AB " P A P B + ^A^B 

Sq-ab - p b <Ja + Pa<5b <8) 
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The T/D 100 then puts the real component into the channel A output 
buffer and the imaginary component into the channel B output 
buffer. 

Probability Density . 

On the T/D 100 the Probability Density of a time -domain in- 
put is ascertained by performing an Amplitude Histogram on the 
input function. This algorithm counts the number of times that 
each of 255 discrete amplitudes occurs at the input. Zero-ampli- 
tude occurrences appear mid-scale of the output display with 
negative-amplitude occurrences increasing to the left of mid-scale 
and positive-amplitude occurrences increasing to the right. 


Average rms (^) 

The average rms of the input signal is calculated from the 
PSD coefficients £s(kf 0 )J by the following equation: 


V 2N 


S(kf 0 ) 


where f 0 » is the base frequency or analysis bandwidth. 
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SECTION 7 


DISCUSSION OF RESULTS 

This discussion will review the major findings of the inves- 
tigation and illustrate each of these findings with appropriate 
examples of the stochastic analysis results. 

The major findings of this study are as follows: 

1. The structural response of the aircraft to buffet 
during moderate-g and high-g maneuvers is very 
complex. Almost all natural vibration modes of 
the aircraft (symmetric and antisymmetric) can 

be excited to a significant intensity level at 
some time during the maneuvers. 

2. As a consequence of 1, above, several different 
types and locations of instrumentation are needed 
to obtain an adequate description of the struc- 
tural responses. 

3. Some of the sensors show pronounced changes in 
the relative modal contributions to the total 
response as penetration (angle of attack) beyond 
buffet onset increases. 

4. The fluctuating shear and bending moment loads 
for the test conditions examined are small in 
terms of design loads except near the wing tip. 
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5. The wing structural response in torsion is much 
larger than anticipated on the basis of previous 
buffet studies , amounting to between 1/3 and 
1/2 of the rms values of wing bending response 
at high angles of attack. 

6. There is some evidence of bending- torsion 
coupling starting at angles of attack between 
9 and 12 degrees at Mach numbers of 0.80 and 
0.87. 

7. The stochastic analysis techniques employed in 
this investigation provide useful and consis- 
tent sets of data despite the short sample 
times and the abruptness of some of the maneu- 
vers. Proper interpretation of the analysis 
results must include consideration of the ranges 

of angle of attack, Mach number and dynamic pres- 
sure that occur during the time duration of each 
record. 

Characteristics of the Structural Responses 

During the course of the study approximately 800 power spectral 
density plots were obtained. These data have two primary uses. 
First they permit identification of the significant modal contri- 
butions through comparison with the ground vibration test data 
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obtained in the F-lll development program. Second they provide 
the data base for future assessment of prediction methods. 

Natural Vibration Modes 

A summary of measured natural vibration inodes of the F-111A 
at 26 degrees leading-edge sweep is presented in Table 8. Most 
of the values were taken from tests previously conducted on Air- 
plane No. 12 which is structurally the same as Airplane No. 13 
(Reference 18) . A few of the measurements were taken on Airplane 
No. 1 and were not repeated on Airplane No. 12 because the parts 
and associated structure had not changed (Reference 19). Data 
are presented for two fuel loading conditions: fuselage empty, 

wing empty and fuselage full, wing empty. Note that the differ- 
ence in fuselage fuel load affects some of the natural frequencies. 

Table 9 presents calculated values of symmetric vibration 
mode frequencies obtained using a math model of the aircraft 
structure for two gross weight conditions which represent specif- 
ic loadings occuring on maneuvers investigated in this study. 

These calculated mode frequencies serve the purpose of clarifying 
some measured responses which were not clearly identified in the 
ground vibration tests because of the close proximity of other 
modes. 
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TABLE 8 - MEASURED F-lliA NATURAL VIBRATION MODES 
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Power Spectra 

Representative examples of power spectra will be presented 
and discussed to illustrate several of the major findings of the 
investigation. All of the plotted PSD's are presented in Appen- 
dix B and tabulated PSD data are presented in Appendix C for 
future comparisons with predictions. 

The power spectra are presented on a consistent semi- log 
format in which normalised values of power spectral density are 
plotted against frequency over the range from 0 to 100 Its (see 
Figure 10(a)). The normalizing factor (scale factor) in each 
case is the sum of the spectral estimates at each frequency from 
1 to 100 Hz. The scale factor is presented above each plot in 
both S.I. units and U.S. customary units (if appropriate). The 
values plotted at 0 and 1 Hz are ficticious and are used only to 
set the scales for the automati - plotting routine that was used 
to generate these presentation-. Any data point which falls on 
the lower bound of the plot at frequencies other than 1 Hz actu- 
ally represents a value below the threshold of the dynamic range 
of the analysis system for the particular analysis. Finally, 
although all the power spectra estimates are plotted from 2 to 
100 Hz, the valid upper limit of frequency for each instrumenta- 
tion item and flight condition is approximately that indicated 
previously in Table 5 because of the flight recorder frequency 
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response characteristics. Values for frequencies much above the 
Units indicated in Table 5 probably are caused by cross-talk 
between recorder channels and should be ignored. 

Near Buffet Onset 

The particular set of PSD data presented in Figure 10 is for 
point 8 of Flight 77, Run S&C-R- The nominal angle of attack is 
7.1 degrees and the length of the sample record is two seconds. 
The aircraft is slightly beyond buffet onset at this condition. 
The wing shear strain-gage data are presented first for each of 
the four wing stations, followed by the wing bending gage and 
torsion gage data. The various accelerometer data complete the 
presentation for this point. Each plot is annotated with the 
vibration node at which each of the major responses occur. 

Examinati.a of the shear data shows that at Wing Station 1 
the major responses are in the first and second wing bending 
modes, the first fuselage bending mode, the first wing torsion 
mode and i combined wing- tail mode. Further outboard on the wing 
some additional modes cam be identified including the second 
antisymmetric wing bending rode and several horizontal tall pitch 
and torsion modes. In addition, the levels of the power spectra 
at the higher frequencies increase to the point that the first 
wing bending and first fuselage vertical bending mode are no 
longer the dominant modes at Wing Station 4. 
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The responses shown by the wing bending gages at each wing 
station track the responses shown by the shear gages at the 
corresponding stations, except that at the inboard two stations 
the first wing bending mode dominates the responses. 

The wing torsion responses are strong at the first wing 
torsion frequency with a major contribution occurring also at the 
first wing bending frequency. There is some response shown at a 
frequency of about 20 Hz which can be attributed to a fuselage 
bending mode. 

Ihe wing tip accelerometers show responses in most of the 
modes already identified from the strain gage sensors plus sever- 
al responses at higher frequencies between 40 and 80 Hz. Note 
that while the PSD p ,ts for the left-hand and right-hand wing 
tip accelerometers are similar they are not identical in detail. 

The fuselage mounted accelerometers show very little re- 
sponse in the first wing tending modes. The two C.G. vertical 
accelerometers show very similar responses and the dominant re- 
sponse is in a band of frequencies from 40 to 45 Hz which in- 
cludes some higher order horizontal tail pitch modes. The other 
significant responses are attributable to first fuselage vertical 
bending, wing-tail, horizontal tail second bending, first 
antisymmetric wing torsion and lower order horizontal tail pitch 
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modes. The dominant responses shown by the pilot station verti- 
cal accelerometer include the wing-tail, horizontal tail second 
bending, first wing torsion, and lower order horizontal tail 
pitch modes. Only small responses occur at first wing bending 
and first fuselage vertical bending modes. The lateral accelero- 
meters at the C.G. and pilot's seat respond at frequencies for 
first and second fuselage lateral bending, antisymmetric wing- 
tail, first wing torsion, second wing antisymmetric bending, and 
horizontal tail pitch modes. Responses are also noted at high 
frequencies between 60 and 80 Hz but these are not readily 
identifiable modes. 

Near Maximum Buffet 

From the data just presented it can be seen that the aircraft 

responses just beyond buffet onset are quite complex and many 

vibration modes are excited. It is observed that a higher angle 

of attack condition where the rms values of response are much 

larger, the relative model contributions change significantly. 
Figure 11 presents the power spectra for point 10 of 

Flight 77, Run S&C-R which corresponds to a nominal angle of 

attack of 12.2 degrees. Examination of the shear data at Wing 

Station 1 reveals a very large response in first wing bending, 

some contributions in the wing second bending mode and new major 

contributions in the first wing torsion and second antisymmetric 

wing bending mode. At Wing Station 2 the contributions due to 
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first wing bending mode and the second antisymmetric wing bending 
mode predominate with significant response occurring in two hori- 
zontal tail pitch modes at higher frequencies. The response in 
the wing- tail and r°cond wing bending modes are now small com- 
pared to the other modes. At Wing Station 3 the first wing bend- 
ing response predominates with some response shown in wing second 
bending and higher order horizontal tail modes. Relatively small 
response is noted at t'.e wing first torsion and the second wing 
antisymmetric bending frequencies. At Wing Station 4 the peak 
power value occurs at the second wing antisymmetric bending fre- 
quency with somewhat smaller responses in first wing bending, 

wing- tail, second wing bending, first wing torsion, and the 
horizontal tail pitch modes. 

The wing bending gage data at the four wing stations show 
responses that are somewhat different from the shear gages. At 
Wing Station 1 the first wing bending mode predominates with 
small responses occurring in first wing torsion, second antisym- 
metric wing bending, and two horizontal tail pitch modes. At 
Wing Station 2, the first wing bending response is very strong 
and a small response is shown for second wing bending. Very 
little response is indicated at higher frequencies. At Wing Sta- 
tion 3, the first wing bending mode and the second antisymmetric 
wing bending mode responses are both large; and significant re- 
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sponses occur at second wing bending, first wing torsion, and the 
horizontal tail pitch modes. At Wing Station 4 the predominant 
mode is the second antisymmetric wing bending mode and smaller 
responses occur at first wing bending, second wing bending, first 
wing torsion, and the horizontal tail pitch mode frequencies. 

A most significant finding of this study is indicated by the 
wing torsion responses. At Wing Stations 1, 2 and 4 a small re- 
sponse is noted in first wing bending* and the major response is 
in a broad band of frequencies varying from about 16 to 35 Hz 
with the major peak corresponding to first wing torsion at Wing 
Stations 1 and 2 and to second antisymnetric wing bending at Wing 
Station 4. This type of response is indicative of a very strong 

bending- tors ion coupling. No torsion response data were available 
at Wing Station 3 for this flight, however, the response is ex- 
pected to be similar in modal content to those at the other wing 
stations on the basis of responses noted for Flight 79, Run 9R. 

The wing tip accelerometer responses for point 10 indicate 
that the left wing and right wing may not be receiving the same 
type of excitation. The left wing tip accelerometer shows a 
broad band of response b- tween 40 and 80 Hz, a fairly large peak 
in a horizontal pitch mode at about 37 Hz and smaller peaks in 
the first wing bending mode and a lower order horizontal tail 
pitch mode. The right-hand wing tip accelerometer response shows 
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peaks in a wing- tail mode, second wing bending and second anti- I 

symmetric wing bending mode in addition to the peak in first wing 
bending and the broad band of responses between 40 and 80 Hz. 

The responses in first wing torsion are not large for either tip ' 

i 

accelerometer which is somewhat surprising in view of the strong 
response indicated at Wing Station 4 on the right-hand wing. 

The two vertical C.G. accelerometers show about the same * 

responses with the major excitation occurring at the second anti- 

f 

symmetric wing bending mode and the horizontal tail pitch modes. 

Some significant response occurs at higher frequencies between 
70 and 100 Hz. The pilot's seat vertical accelerometer shows a 
moderate response in a band including first wing bending and 

{ 

f 

first fuselage vertical bending, and major responses in second 
wing bending, first wing torsion, second antisymmetric wing bend- 
ing, and the horizontal tail pitch modes. The largest peak is at 
the second antisymmetrical wing bending frequency. 

* 

The lateral accelerometer at the pilot's seat has majo- 
responses at the second fuselage lateral bending, second anti- 
symmetric wing bending, and a high order horizontal tail mode 
with only moderate responses in the first fuselage lateral bend- 
ing mode and a wing- tail mode. The responses of the C.G. lateral 
accelerometer are similar in modal content to those at the pilot's 
seat except the wing-tail and second fuselage lateral bending 
mod 3 responses are much smaller. 1 
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Interpretations 


The right wing exhibits large responses in several inodes and 
in particular undergoes a strong bending- tors ion coupling. The 
C.G. and pilot seat accelerometers respond primarily at frequen- 
cies above 16 Hz despite the fact that large wing responses are 
noted at lower frequencies. A significant amount of the wing 
excitation occurs in antisymnetric modes despite the fact that 
the maneuver is nominally a symmetric maneuver. Further inves- 
tigation brought out the fact that points 9 and 10 each include 
a brief roll disturbance and recovery. Figure 12 shows the 
roll-rate time history of the maneuver under study. After initial 
entry into the wind-up turn the roll- rate decays in a slightly 
oscillatory manner and just after 18 seconds a small spoiler in- 
put drives the roll rate positive. The rate decays rapidly to 
near zero and then just after 22 seconds another small spoiler 
input again drives the roll rate positive. Referring back to the 
other time histories in Figure 5(b)it can be seen that each spoiler 
input is preceded by a change in character of the torsion re- 
sponses. This fact is particularly evident just before 22 sec- 
onds. In the sideslip angle trace in Figure 5(e)it can be seen that 
some oscillation occurs. All of these indications point to the 
occurrence of mild wing rock which possibly induces the antisym- 
metric responses noted in Figure 11. 
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Figure 12. ROLL RATF TIME HISTORY, FLIGHT 77, RUN S&C-R 



It is appropriate now to review information learned so far 
First, the structural responses during maneuver well into buffet 
are complex. Many of the natural vibration modes of the aircraft 
are excited and develop significant responses. Second, several 
different types of sensors and sensor locations are required to 
obtain an adequate description of the aircraft structural re- 
sponses. Third, the modal content of the responses can change 
significantly as the penetration beyond buffet onset increases. 

As a consequence of these three findings it appears unlikely 
that a simplified prediction method will produce accurate 
estimates of buffet loads and accelerations, but there is a need 
to assess how severe the buffet loads and accelerations are before 
making any judgment on what degree of accuracy may be needed. 

Magnitudes of the Structural Responses 
The complexity of the modal responses makes it difficult to 
comprehend the variations in magnitude of the structural responses 
by attempts to compare mode by mode. Therefore, the concept of 
root -mean- square (rms) values is used for this purpose. By summing 
the magnitudes of the power spectra over a range of frequencies 
and taking the square root of the sum, a measure of the total 
response for each item of instrumentation is achieved. 

In the following discussion the rms values are evaluated over 
the frequency ranges from 1 to 50 Hertz or from 1 to the limit of 
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recorder response if less then 50 Hertz for a particular item. 

If nns values over a different frequency range are desired they 
can be calculated using the tabulated PSD data presented in 
NASA CR 152111. 

One purpose of the rms analysis was to investigate effects of 
Mach number and altitude on the magnitudes of the responses. To 
do so, first the results of the wind up turn at M«0.70 are examin- 
ed, followed with comparisons of the data for the three M«0.80 
maneuvers, and lastly, with comparisons of the three maneuvers 
at M-0.87. 

M*0.70 Results 

Figures 13 through 15 present the wing strain gage responses 
at each of the four instrumented wing stations. The rms values 

i 

are plotted at the nominal angles of attack previously assigned 
in Table 7 (remember, however, that each data point represents 
data taken over a discrete interval of time and therefore a range 
of angles of attack). The curves shown in Figures 13 through 
15 show a remarkable consistency in the variations. In every 
case the trend with angle of attack is similar. In particular 
an inflection occurs between 11 and 12 degrees nominal angle of 
attack. Referring back to Figure 6(a) it can be seen that an 
inflection also occurs in the variation of angle of attack with 
time. The implication is that a wing panel leading-edge stall 
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has occurred. Figures 13 and 14 show that the rms values of 
shear and bending moment decrease rather abruptly beyond 11 degrees 
and then rise again somewhat. Figure 15 shows that the torsion 
values tend to remain nearly constant above 11 degrees. Inspection 
of the PSD data for this maneuver (NASA CR- 152110) revealed 
that between 11 and 12 degrees there is a shift in bending moment 
response from the first wing bending mode to higher frequency 
bending modes, particularly outboard on the wing. Finally the 
time histories for the right-hand wing spoilers presented in the 
Appendix show that a small and short duration spoiler deflection 
occurs just after the time sample associated with the data point 
represented at 11 degrees angle of attack. All this evidence 
points to a mild wing rock occurring during development of the 
leading-edge stall. 

Figure 16 presents the accelerometer responses obtained 
during the wind up turn maneuver of Flight 48, Run 6. Figures 
16(a) throjgh 16(d) show the vertical accelerometer responses at 
the right-hand wing tip, the center of gravity and the pilot's 
seat. The similarity of the trends in variations of the re- 
sponses with angle of attack is again noted. All the vertical 
accelerometers show an inflection occurring between 11 and 12 
degrees just as the wing strain gage responses did. The struc- 
tural characteristics of the aircraft are such that a significant 
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reduction in rtns magnitudes of the response occurs between the 
wing tip and the center of gravity and a further reduction occurs 
between the center of gravity and the pilot's seat. 

Figures 16(e) through 16(h) present the vertical accelero- 
meter responses for both the right and left hand wing tip and 
the lateral accelerometer responses at the pilot's seat and the 
center of gravity. Comparing Figures 16(e) and 16(f) it is seen 
that the responses at the two wing tips are similar at the lower 
angles of attack and then differ somewhat at the higher angles 
of attack. The differences may be attributed to a slightly 
different progression of flow separation on the two wings which 
is accentuated by the spoiler deflection on the right hand wing 
between the data points shown for 11 and 12 degrees angle of 
attack. 

Figures 16(g) and 16(h) show that the lateral accelerations 
are small at the pilot's seat and particularly so at the center 
of gravity. 

M ■ 0.80 Results 

Figures 17, 18, and 19 present the wing strain gage re- 
sponses and Figure 20 the accelerometer responses for the three 
maneuvers performed at approximately M - 0.80 but at different 
altitudes. Four sets of data are shown because both one-second 
and two-second samples wer? analyzed for the wind up turn maneu- 
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ver of Flight 77, Run S&C-R. These are shown on the plots as the 
circles and squares, respectively. These results will be discussed 
first before considering the responses for the other maneuvers. 
While the one-second and two-second data are generally similar 
in magnitude, some differences do exist particularly in the high 
angle of attack region at about 12 degrees for the shear and 
bending moments. The fact that a mild wing rock occurs during 
this maneuver has previously been discussed under the power 
spectral density analysis. It is apparent that the two-second 
samples tend to smooth the variation of response with angle of 
attack for Wing Stations 1, 2 , and 3, but there is little dif- 
ference at Wing Station 4. The torsional responses tend to be 
more nearly the same for the one- second and two-second sample 
data at each angle of attack. 

The general trends of the variations in response with angle 
of attack are different from the trends at H * 0.70 in two ways. 
First buffet onset occurs at lower angle of attack (about S de- 
grees) which is expected for this particular wing leading-edge 
sweep. The initial flow separation is shock- induced at this 
condition. Second, the increase in the responses with angle of 
attack is generally more gradual even though the maximum measured 
rms values are much larger than those at M ■ 0.70. Comparison 
of the shear and bending moment response magnitude variations 
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with wing station indicate that the flow separation is primarily 
outboard of Wing Station 2, since there is little change in rms 
values of shear between Stations 1 and 2 but a large change in 
rms values of bending moment between those stations. The torsion 
responses at Wing Stations 1, 2, and 4 all have maximums at about 
12 degrees angle of attack and then decrease, whereas at M * 0.70 
the torsion data indicated a nearly constant response level at 
high angle of attack. Torsion data are not available at Wing 
Station 3 for Flight 77. 

The vertical accelerometer responses shown in Figures 20(a) 
through 20(d) show little difference in the right-hand wing tip 
accelerations for the one- and two-second samples, fairly signif- 
icant differences for the center of gravity measurements and 
smaller differences for the pilot's seat measurements. 

The wing tip responses are three times larger at M ■ 0.80 
than they were at M • 0.70, whereas the center of gravity re- 
sponses and the pilot seat responses are about twice larger. 

The variations with angle of attack are similar to each other 
with peak magnitudes occurring at about 12 degrees. 

Figures 20(e) and 20(f) show the accelerations a. the right- 
hand and left-hand wing tips in a side-by-side presentation to 
allow comparison by inspection. Nbte that the variations with 
angle of attack and the peak magnitudes are generally similar for 
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both wing tips for Flight 77, Run S&C-R, but the left-hand wing 
shows somewhat more response than the right-hand wing at the 
highest angle of attack. / 

Figures 20(g) and 20(h) present the lateral acceleration 
responses at the pilot's seat and the center of gravity, respec- 
tively. For this maneuver the responses show similar trends with 
angle of attack and the magnitudes are nearly the same at the two 
locations. The lateral accelerations also peak at about 12 de- 
grees angle of attack, and are about one- fourth of the value of 
the vertical acceleration at the center of gravity and one- third 
the values of vertical accelerations at the pilot's seat. The 
peak magnitudes of the lateral accelerations are about 3 times 
the values at M ■“ 0.70. 

Summarizing the differences between responses obtained dur- 
ing the wind up turn maneuvers at M=0.70 and M"0.80 it was found 
that the differences in flow separation conditions leads to buf- 
fet onset at lower angle of attack at M ■ 0.80, and more gradual 
increases in the responses with angle of attack. The minimum 
measured responses, however, are much larger than those measured 
at M » 0.70 and cannot be explained solely on the basis of the 
increase in dynamic pressure. 

Next, the effects of altitude (or dynamic pressure) on the 
responses are examined. Figures 17 through 19 also present 



wing strain gage data and Figure 20 the accelerations for two 
pull up maneuvers performed at approximately M ■ 0.8 — one at 
about 3800 m and the other at about 1500 m. Data for these 
maneuvers are shown as triangles and diamonds, respectively. In 
general the responses for the lower altitude maneuvers are higher 
for a given angle of attack. A cursory examination of the bend- 
ing moment and torsion responses indicates the magnitudes are 
approximately related between the two pull up maneuvers by the 
dynamic pressure values. Since the power spectra indicate the 
bending moment responses for the two inboard stations are pre- 
dominantly in the first wing bending mode, one can apply the 
simple linear elastic system analysis described in Reference 20. 
That analysis indicates that damping is (1) primarily structural 
if the response relationship is a direct function of dynamic 
pressure or (2) primarily aerodynamic if the relationship is a 
square root function of dynamic pressure. The measurements dis- 
cussed above would indicate that structural damping predominates 
in the first wing bending mode for the two high-g pull up maneu- 
vers at M ■ 0.80. 

It is obvious that the response data for the wind up turn 
do not belong to the same family of curves as the data for the 
two pull up maneuvers. There are several pertinent factors which 
can alter the response such that the expected relationship did not 
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apply. First the difference in pitch rate between the wind-up 
turn and the pullups is a possible cause. Second, the gross 
weight for the wind-up turn was much lower than those for the 
pullups. Third is the fact that the wing has a different 
aeroelastic shape (twist distribution) for each maneuver altitude 
and the progression of flow separation can be significantly diff- 
erent. This aspect of the problem will be explored further as 
we examine the higher Mach number results. 

M — 0.87 Results 

Figures 21 through 24 present the wing strain gage and 
accelerometer responses for the three maneuvers performed at Mach 
numbers of about 0.87. Because these maneuvers were all rather 
abrupt, the rms magnitudes of response were expected to be similar 
in trend with angle of attack and to correlate with altitude. How- 
ever, that did not happen. The results of the roller coaster 
maneuver performed at high altitude will be examined first, follow- 
ed by further examination of the effect of Mach number on the var- 
iations of response with angle of attack. 

The maximum magnitudes of shear and bending moment response 
are about the same as those measured at M"0.80 but the values 
for torsion response are somewhat lover. There is an indication 
that a peak might occur in some of the responses just above the 
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Figure 21. ROOT MEAN SQUARE VALUES OF WING SHEAR RESPONSES 
FOR NOMINAL MACH NUMBER - 0.87 
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Figure 22. ROOT MEAN SQUARE VALUES OF WING BENDING MOMENT RESPONSES 
FOR NOMINAL MACH NUMBER -0.87 
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Figure 24. ROOT MEAN SQUARE VALUES OF ACCELEROMETER RESPONSES 
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Figure 24. CONCLUDED 




last data point or at about 15 degrees angle of attack. The 
accelerations are in general lower than at M * 0.80 and do not 
peak in the angle of attack range measured. 

Turning now to the data for the pullup performed at about 
3700 m which are represented by square symbols in Figures 21 
through 24, it can he seen that, in general, the variations in response 
are similar to those for the roller coaster maneuver. The higher 
initial levels of response are caused by the fact that the wing 
was in buffet even at trim conditions for 1-g flight. It is 
interesting that the change in slope of the responses still 
occurs above 9.3 degrees angle of attack. At the higher angles 
of attack the magnitudes of response for the shear and bending 
moments are slightly smaller than those for the roller coaster 
maneuver despite the fact that the dynamic pressure is much 
larger. The fuselage vertical accelerations follow a similar 
trend but the wing tip accelerometer responses are slightly higher 
than those for the roller coaster maneuver. 

Data for the pullup maneuver performed at 1500 m are pre- 
sented as triangles in Figures 21 through 24. Because this 
maneuver was so abrupt only three points are shown, one at 1-g 
trim conditions and two at high angle of attack. The wing strain 
gage responses at high angle of attack are larger than can be 
accounted for by a dynamic pressure effect. The probable cause 
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is the fact that the Mach number is about 0.04 lower than the 
other maneuvers at this angle of attack. The accelerometer re- 
sponses are only slightly higher than those for the other pullup 
maneuver . 

When the responses for the three maneuvers are compared, it 
is again obvious that they do not correlate with altitude. In this 
comparison, two of the three factors mentioned in the discussion 
of the M«0.80 comparisons have been minimized. The pitch rates 
are more nearly the same for the present comparison and the gross 
weights are also nearly the same. As a consequence we must 
attribute the anomalous trends with altitude primarily to differ- 
ences in aeroelastic deformation of the wing which can cause 
significantly different shock patterns and thus differences in 
the progression of flow separation with angle of attack. 

Summary Analyses 

Before moving on to a discussion of the significance of the 
measured responses with respect to design loads, a brief review of 
what has been learned about the effects of Mach number is in order. 

Figure 25 presents comparisons of the shear, bending moment and 
torsion responses at Wing Station 1 and the pilot's seat verti- 
cal accelerometer responses as functions of angle of attack 
for each of the three high altitude maneuvers. Three features 
of these comparisons are of particular interest. 
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Figure 25. TYPICAL EFFECTS OF MACH NUMBER ON WING 
AND PILOT'S SEAT RESPONSES 




First is the different response variations with angle of 
attack as Mach number increases. At M * 0.70 the variation is 
typical of that produced by leading-edge flow separation, which 
starts outboard on the wing and then progresses rapidly across 
the wing as angle of attack is increased. The peak response 
apparently occurs when the progression is completed. At M = 0.80 
the variation is representative of shock- induced separation under 
conditions where the outboard shock is located very near the 
crest of the airfoil and moves ahead of the crest as angle of 
attack is increased. At M = 0.87 the variation is representative 
of a shock induced separation in which the outboard shock is well 
aft of the crest at low angle of attack and then moves forward 
slowly with angle of attack until the shock passes forward of the 
crest. 

Second is the fact that the torsional responses tend to peak 
and the angle of attack for the peak increases as Mach number in- 
creases. This behavior may be a manifestation of the breakdown 
of the flow from a shock- induced separation to a leading- edge 
separation at the higher Mach numbers as described by Pearcey 
in Reference 21. 

Third is the fact that the pilot’s seat vertical accelero- 
meter response tends to track the wing torsion response more so 
than it does wing shear or bending responses both in terms of 
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curve shape and relative magnitudes. This finding could lead to 
a whole new approach to predicting pilot's seat response if the 
same trends continue for other wing sweep angles. In any event, 
it indicates that scaling of the first wing torsion frequency, in 
addition to the first wing bending frequency, may lead to more 
accurate predictions of full-scale responses from small-scale 
wind tunnel model measurements. 

Normalized Accelerometer Responses 

In order to gain a perspective of the relative magnitudes of 
the accelerometer responses during buffet, the maximum measured 
responses were normalized by the maximum maneuver load factor for 
each maneuver. The normalized responses are presented in Figure 
26. The various symbols designate nominal values of altitude and 
the Mach number values represent the average Mach number for each 
particular data sample. The analysis shows three pertinent 
results. 

First, the relative responses at high altitude are larger 
than those at lower altitudes. This result was expected since 
the aircraft must operate at higher angle of attack for a given 
load factor and therefore penetrates farther above buffet onset 
at high altitude than it does at low altitude. 
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Second, the maximum relative response for the wing tip and 
pilot's seat vertical accelerometers increases with Mach number, 
reaches a peak between 0.80 and 0.85 Mach number and then de> 
creases again whereas the center of gravity response continues to 
increase at high altitude. As a consequence, the center of 
gravity response is probably not the best measurement to use as 
an indicator of buffet intensity effects. 

Third, the structural arrangement of the aircraft is such 
that response levels attenuate considerably between the wing tips 
and center of gravity and also attenuate between the center of 
gravity and the pilot’s seat. A major factor contributing to the 
attenuation is the fact that very little response is excited at 
those two fuselage locations at the low frequencies associated 
with the first symmetric and first antisymmetric wing bending 
modes with the wing sweep set at 26 degrees. 

Normalized Wing Buffet Loads 

One of the key questions that was to be answered by this 
investigation was just how significant are the buffet loads 
encountered in high-g maneuvers. It was possible to obtain an 
answer to this question for the wing shear and wing bending loads 
because M*,88, ^ LE -26 was one of the anticipated limit loading 
conditions for the structural design of the aircraft and detailed 



analyses of the "steady" (mean) loads produced during the Loads 

Demonstration Program for the F-111A aircraft. The wing torsional 

moment at this same condition is not a critical loading for the 

F-111A aircraft. The maximum rms values of shear and bending 

moment measured at each wing station during the three maneuvers 

performed at about M=0.86 have been normalized on the basis of 

the design load factor "steady" airloads for the specific dynamic 
pressure values associated with the measured buffet loads. The 

"steady" airloads were obtained from an analysis presented in 

Reference 22 in which values of shear and bending moment per unit 

dynamic pressure were plotted against dynamic pressure for ^“26°, 

M*0.88. 

The results of this analysis are presented in Figure 27 
which shows the spanwise variations of normalized maximum shear 
and bending moment buffet loads. The symbols designate the 
different maneuvers. The plots show that there is an important 
spanwise variation of buffet loading and that the highest rela- 
tive values occur at the most outboard measuring station. The 
effect of altitude is essentially similar to that for the accel- 
erometer responses. At high altitude the dynamic response is a 
larger percentage of the "steady" airload than at lower altitudes. 
While the relative buffet loads are small they are not insignif- 
icant. 
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SECTION 8 


CONCLUDING REMARKS 

This investigation has provided some valuable new information 
about aircraft structural response to aerodynamic buffet at sub- 
sonic and transonic speeds. In particular, it was found that the 
structural response to buffet during moderate- to-high-g maneuvers 
is very complex. Nearly all the natural vibration inodes of the 
aircraft, both symmetric and antisymmetric, are excited during a 
maneuver. As a consequence, an extensive array of instrumentation 
is needed to obtain an adequate description of the structural re- 
sponses. As the penetration (angle of attack) beyond buffet on- 
set increases, significant changes in the modal contributions to 
the total roct-mean-square responses can occur. Therefore, no 
single vibration mode can be selected as an indicator of the 
variation of buffet intensity during a maneuver. 

In general, the fluctuating shear and bending moments on the 
wing are small except near the wing tip where the rms values 
approach 15 to 20 percent of the design "steady" loads. The rms 
values of torsion response are higher than expected on the basis 
of previous buffet studies and amount to between 1/3 and 1/2 of 
the corresponding wing-bending response values. Some evidence 
of strong wing bending- torsion coupling was obtained at high 
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angle of attack at transonic flow conditions. 

The effects of Mach number on the variations of wing responses 
with angle of attack were about as expected from previous buffet 
studies. At M * 0.70 the progression of response intensity with 
angle of attack is quite rapid and reaches a peak about four de- 
grees after onset and then decreases. At M * 0.80 onset occurs 
at a lower angle of attack than at 0.70. The progression is 
somewhat less with angle of attack but the maximum measured 
values are significantly larger and occur at the highest angle 
of attack. At M « 0.87 onset occurs at a still lower angle of 
attack than at 0.80. The initial rate of increase of response 
is low for several degrees and is followed by a larger rate of 
increase such that the maximum measured values are about the same 
as those measured at M ■ 0.80. 

Attempts to correlate the responses at the various alti- 
tudes as functions of dynamic pressure generally failed. 

There is an indication that damping in the first wing bending 
mode is primarily structural at low altitude but not at high 
altitude where the "steady” air loads are much smaller in ab- 
solute magnitude. On the basis of the limited data available 
it was concluded that the differences in "steady" aeroelarftlc 
deformation (twist) of the wing alter the progressior. of flow 
separation as angle of attack is increased at the different 
altitudes. 
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J. G. Jones has shown in Reference 23 that the ’'apparent 
damping" in the first wing bending mode is highly nonlinear with 
angle of attack for a small fighter aircraft. It is possible 
then that the variations in the aerodynamic damping with angle of 
attack are different for the F-lll wing at different altitudes 
due to the differences in twist distribution. That situation 
could account for the inability to correlate the effects of 
altitude with functions of dynamic pressure. 

It is now appropriate to review what impact the results of 
this investigation might have on future flight programs to inves- 
tigate buffet characteristics and also on the development of 
methods to predict structural responses to aerodynamic buffet. 

First, with respect to future flight programs, it is clear 
that a significantly larger array of instrumentation is needed to 
adequately assess aircraft structural response than has been used 
on most previous flight programs. In addition the sensitivities 
of strain-gage-bridges used to measure buffec loads should be 
significantly larger than those used to measure total aircraft 
loads in order to obtain an optimum signal- to-noise ratio for 
processing the dynamic data. In-flight measurement of the 
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fluctuating pressures and the aeroelastic deformations of the 
wing could considerably enhance the value of flight measurements 
of buffet characteristics that are to be compared with predictions. 
The twist distributions are of primary interest although chordwise 
deflections at several span stations would be of significant 
value if fluctuating pressure measurements are obtained in addi- 
tion to loads and accelerations. 

All measurements that are pertinent to a buffet analysis 
should be recorded on FM tape recorders and subchannel recording 
frequency limits should be no lower than the second wing torsion 
frequency. If leading and/or trailing edge flaps are a part of 
the configuration, frequencies associated with flap modes should 
also be within the recorder frequency limits. Finally, it is reco- 
mmended that very gradual windup turn maneuvers be used to acquire 
buffet measurements so that sample times will be adequate for 
statistical analyses. 

Second, with respect to development of prediction methods, 
one must consider the objective of the prediction. If it is 
desired to obtain detailed estimates of the aircraft vibrational 
environment due to buffet, i.e., loads, accelerations, and dis- 
placements, then the prediction method should probably include 
considerations of all aircraft symmetric and antisymmetric vibra- 
tion modes up to and including that associated with second wing 
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torsion. If the objective is to obtain comparative rms accelera- 
tions at the pilot’s seat at different flight conditions, it 
would probably be sufficient to consider symmetric and antisym- 
metric modes up to and including first wing torsion and symmetric 
and antisymmetric wing-fuselage-horizontal tail modes in the 
range of frequencies covered by the wing modes. The horizontal 
tail modes appear to be particularly important with respect to 
damping of the fuselage motions. 
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APPENDIX 


TIME HISTORIES OF INSTRUMENTATION OUTPUT 

In this Appendix time histories of the Instrumentation output 
are displayed in analog form for the following maneuvers: 


Figure 

Flight 

Run 

Maneuver 

A-l 

48 

6 

Windup Turn 

A- 2 

78 

5 

Pullup 

A- 3 

79 

9R 

Pullup 

A* 4 

60 

10 

Roller Coaster 

A- 5 

78 

4 

Pullup 

A- 6 

70 

2 

Pullup 


The time histories for the Flight 77, Run S&C-R windup turn maneu- 
ver are presented in the main body of the report. 

Each set consists of several records and each record contains 
the output from several instrumentation channels. The top channel 
on each record shows the output for the same item - AW002, the R/H 
wing tip accelerometer (except Figure A- 3 which has DH001C, angle 
of attack). The particular groupings of items presented are those 
used in the study and seme items appear more than once. The flight 
time (actually frame number) which is the correlating item appears 
at the lower boundary of each record. Each second is indicated 
by a tick in the time trace. 
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Figure A-2. TIME HISTORIES OF INSTRUMENTATION OUTPUT 
FOR PULLUP, FLIGHT 78, RUN 5 



Figure A-2. Continued 
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Figure A-4. Continued 


167 















I 


original page is 

POOR QC/-.UTY 



u V U.- AaA 





!li 


*> *it»* lf.7f.alir 


<N) (IRS) 


*1 M»m » * _'»« !«• 


♦4 ’A*in 1 1.0/* 



* lo *«. i7*ro 


f N) flfts) 


(mr> (wr> 


(i*o (wr> 


> t ♦ l i w S 


t*Au /kao'** *i ) 


-4- — — ; 

)€ 57 5» 5* 40 

Tfl« — f VOMDS 


Figure A -5. Continued 


* 'S 



ifc« i AwJLli jLb a i 



» 1*4*4 ** \fc \? *« 

ni* — sfoarvc 


(d) 


Figure A-5. Continued 
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Figure A-5. Concluded 




t ut r+i mu m 


iSS 


%W*>I 

n-T 

\*«j 


: n n n 

: i t nr 


>41 . 

U v 
» f 


4 •»»!*- 
vn'i 



s I 
tmi r . 


TW - »tU>NK 


I '> hh 
♦ I Na Jl*** 


I l»» 111 

\ 

I1H» ill 


2 uniu 
♦4 


f» *41 


I 1 *■ 

,M N 


n 

» I H.'t IU 


I M/k * • 
•4 4i*;. > 


-4 41*1*1 

♦ i a ;*!!! 1 


n*i*i 


I 4?* It* 
♦ / 


(a) 


Figure A-6. TIME HISTORIES OF INSTRUMENTATION OUTPUT 
FOR PULLUP, FLIGHT 70, RUN 2 
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Figure A-6. Continued 
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Figure A-6. Concluded 
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